Detection of transplanckian effects in the cosmic microwave background 
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Quantum gravity effects are expected modify the primordial density fluctuations produced during 
inflation and leave their imprint on the cosmic microwave background observed today. We present 
a new analysis discussing whether these effects are detectable, considering both currently available 
data and simulated results from an optimal CMB experiment. We find that the WMAP data 
show no evidence for the particular signature considered in this work, but give an upper bound 
on the parameters of the model. However, a hypothetical experiment shows that with proper 
data, the transplanckian effects should be detectable through alternate sampling methods. This 
fuzzy conclusion is a result of the nature of the oscillations, since they give rise to a likelihood 
hyper-surface riddled with local maxima. A simple Bayesian analysis shows no significant evidence 
for the simulated data to prefer a transplanckian model. Conventional Markov Chain Monte Carlo 
(MCMC) methods are not suitable for exploring this complicated landscape, but alternative methods 
are required to solve the problem. This however requires extremely high-precision data. 



INTRODUCTION 



Quantum fluctuations in the scalar field responsible 
for driving inflation can give rise to a power spectrum of 
primordial density perturbations consistent with what is 
required to seed the large-scale structure (LSS) and the 
temperature anisotropics in the cosmic microwave back- 
ground (CMB) radiation. This fact is probably one of the 
main reasons for inflation having become a part of the 
concordance model of cosmology today. The quantum 
fluctuations are imprinted during the early inflationary 
epoch, and become classical fluctuations in the gravita- 
tional potential as the fluctuations leave the causal hori- 
zon. The standard way of calculating the power spec- 
trum of the density fluctuations make use of the fact that 
space-time becomes Minkowskian in the distant past. 
However, it has been argued 0, 0, 0, 0, 0, 0, 0, 0, 0| that 
since a given length scale observable in the universe to- 
day shrinks and becomes smaller than the Planck length 
if followed sufficiently far back in time, effects of quan- 
tum gravity will at some point play a role in setting up 
the perturbations, and may potentially leave an observ- 
able signature. The primordial density fluctuations in the 
early universe may have small oscillations superimposed 
as a result of quantum gravitational effects. These effects 
are due to non-negligible curvature imposed by high en- 
ergy densities in dc Sitter space during the early stages of 
inflation. Several papers have suggested a generic shape 
of these effects [!, 0, 0] in the form of small oscillations 
superimposed on the standard, nearly scale-invariant pri- 
mordial power spectrum of density fluctuations. A cru- 
cial questions is, of course, whether such a signature is 
observable in the universe today. This question has been 
investigated in several papers [1, 0, 0, 0, 0, 0, ■ The 
conclusion in (Io| was pessimistic: the so-called trans- 



planckian effects were found to be unobservable in prac- 
tice. However, this conclusion has been called into ques- 
tion in 0] . One of the problems pointed out in that work 
is that the oscillations require more numerical accuracy 
than what is commonly employed in codes like CMB- 
FAST [ll| and CAMB 0|. 

In this paper, we revisit the question of the observabil- 
ity of transplanckian effects in the CMB, using a more 
accurate numerical treatment than in [Iol | . Both present 
data and a hypothetical optimal CMB experiment are 
considered. Furthermore, we mention alternative sam- 
pling algorithms to overcome the numerical problems en- 
countered in exploring what turns out to be a particularly 
nasty likelihood hyper-surface. Section II gives a brief de- 
scription of the particular model we consider and the data 
sets we use. In section III we discuss some technical de- 
tails in the numerical computations undertaken. Section 
IV then contains results obtained with currently available 
data sets, and in section V we investigate the possibilities 
with an ideal CMB data set. The exact likelihood func- 
tion is described in section VI, while a Bayesian evidence 
analysis is performed in section VII. A discussion of the 
results and a suggestion for an improved approach to the 
search for transplanckian effects can be found in section 
VIII, and we conclude in section IX. 



II. BACKGROUND 

Neglecting the perturbations in the metric during in- 
flation, it is viable to assume that the inflaton scalar field 
can be expressed as: 

0(x,t) - 0oW + gM (1) 

classical evolution fluctuations 

The fluctuations in a scalar field 5<f> in a Friedmann- 
Lemaitre-Robcrtson- Walker background satisfy the 
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equation of motion (EOM): 
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+ 3HS(/) - -(V 2 )d> + m z 5<t> = 



(2) 
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where the dot denotes the derivative with respect to time 
and to 2 oc d 2 V/d(f) 2 , where V(cf>) is the inflaton potential. 
Performing a Fourier expansion of the fluctuations: 
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where N is a normalization constant, equation (|2|) be- 
comes: 



linear [Tj] in H/A. In [18|], it is argued that equation ((SJ) 
represents a generic effect, regardless of the details in the 
(undetermined) transplanckian physics. In equation 
([5]) is further parametrized as 



P(e,Z,k)=P {k)[l-t 
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(9) 



with variables explained in table |U Notice how the stan- 



il k + 3Hu k + ( — + m 2 )u k = 0. 



(4) 



Assuming dc Sitter space, this equation can be rewritten 
using a change of variable (uu = ^p/o) and switching to 
conformal time (drj = dt/a). This removes the friction 
term, and the final EOM becomes: 



^pk = 



(5) 



where ' denotes derivative with respect to conformal time 
n. A detailed derivation of this result can be obtained in 
jl3j |. The general solution for ip k is: 
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where A k and B k are Bogoliubov coefficients satisfy- 
ing \Ak\ 2 — \Bk\ 2 = 1, 77 is the conformal time and k 
is the wave mode. The usual choice of vacuum is the 
Bunch-Davies vacuum, which corresponds to A^ = 1 and 
Bk = 0. Space-time becomes Minkowskian in the infi- 
nite past (77 — ► —00), where it is feasible to impose initial 
conditions on the field equations. 

In [H, different possibilities for suitable vacua are dis- 
cussed before focusing on a specific choice, the adiabatic 
vacuum, for which afc(f?o)|0, 770) = 0. This vacuum (for a 
given ?7o) is a typical representative of vacua, and corre- 
sponds to a minimum uncertainty between field and its 
conjugate momentum [Til ]. When imposing initial condi- 
tions, Bogoliubov-coefficients are constrained by the re- 
lation 01: 
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The slow roll parameter 
The transplanckian scale 
The reduced Planck mass 
A parametrization of H/A 
The largest measurable scale 
The standard power spectrum 



TABLE I: Explaining equation © 



dard inflationary power spectrum Pq is restored as the 
cutoff £ oc i — > 0. It is also interesting to note that for 
low cutoff values £ oc -r- — > 00, the amplitude of the mod- 
ulation converges to 2 (as lim^^oo x sin(2/x) = 2). The 
amplitude of the oscillations never reaches values greater 
than 2, and only severely impacts Pq for large values of 
£. Physically, if the cutoff £ oc i was enforced at low 
values (large £), it should easily be detectable in the the 
cosmic microwave background, as will be presented in the 
following section. This is, however, not the case, and we 
know that the transplanckian physics leave imprints at 
£ <C 1, corresponding to a cutoff at high energies. The 
amplitude of the oscillations is therefore expected to be 
very small. 

When evolving the modulated primordial power spec- 
trum through CAMB, the transplanckian imprints for 
varying £ in the full power spectrum can be seen in figure 
[U The oscillations singled out for varying £ are presented 
in figure [2] 
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Note that the Bunch-Davies vacuum is restored in the 
infinite past: When 770 ~~ ► —00, Ak = 1 and Bk = 0. A 
discussion of the implications of the choice of vacuum can 
be found in [IIGJ]. 

When choosing this zeroth order adiabatic vacuum [H , 
the primordial curvature spectrum can be expressed as 



{H\2{H\2{ H ,2A 



(8) 



where A is the (Planck) energy cutoff scale and H the 
Hubble parameter. Note that the size of the corrections is 




FIG. 1: The full power spectrum Ce for varying £. 
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FIG. 2: The transplanckian modulations ACe for varying f. 



A. A brief summary of recent results 

Reference Q concludes that effects of transplanckian 
physics are possibly within the reach of cosmological ob- 
servations. Equation ((9|) describes a a generic expression 
for how transplanckian effects would modify the primor- 
dial power spectrum. As seen in the previous section, 
these oscillations are caused by a nontrivial vacuum for 
the infiaton field 0, E| ■ As the oscillations are expected 
to contribute to the energy density, this could change the 
way the universe expands, and in a worst-case scenario, 
the inflationary phase could be destroyed. The authors 
of [I| investigate this possibility, and conclude that the 
the back reaction is under control and fully consistent 
with inflation, with a slow roll found to be completely 
dominated by the vacuum energy given the parameters 
suggested in [f|. 



1. WMAP data and transplanckian effects 

A previous data analysis [3 concluded that no 
significant signals from transplanckian effects are visible 
in the CMB. Another analysis [|| claimed that there 
are some weak hints in the current data, and these 
indications have become slightly stronger with the 
WMAP3 data compared to earlier claims by the same 
authors [1, d, 0, 0] ■ The parameters implied by the data 
suggests oscillations in amplitude that are periodic in 
the logarithm of the scale of the CMB fluctuations, just 
as predicted from transplanckian physics. 

In Q , the authors discuss the so-called cosmic variance 
outliers, i.e. points which lie outside the ler cosmic 
variance error. These outliers are considered interesting 
as the probability of their presence is very small [2fj| . 
The authors mention that it has been envisaged that 
the outliers could be a signature of new physics, even 
though the cosmic variance could be responsible for 
their presence. The conclusion of [j| is that there exist 
statistical justification for a presence of oscillations in 



the power spectrum. It should however be noted that 
the authors of [|| treat the amplitude of the oscillations 
as a free parameter, which ensures the fitting of the cos- 
mic variance outliers with the transplanckian oscillations. 

Another paper [To| was less optimistic, and concluded 
that it is unlikely that a transplanckian signature of this 
type can be detected in CMB and large-scale structure 
data. However, left some technical issues open, and 
we find it worthwhile to revisit this problem. We show 
that the current WMAP3 data gives weak constraints on 
the parameters of the transplanckian model using con- 
ventional MCMC methods. 



2. Simulated data and transplanckian effects 

The conclusion of was that CMB and LSS data 
are in principle sensitive to transplanckian modulations 
in the primordial power spectrum, but that is practically 
impossible to make a positive detection even in future 
high-precision data. This has to do with the nature of 
the oscillations, as the value of the likelihood function 
is extremely sensitive to e and £. But what "extremely 
sensitive" means was not explained in . We will show 
explicit examples of how this sensitive likelihood function 
contemplates trouble, and why this renders the underly- 
ing MCMC method in CosmoMC of little use in estimat- 
ing the transplanckian parameters. 

III. NUMERICAL DETAILS 

CAMB [l2j was modified to include the primordial os- 
cillations described in equation ©, adding the param- 
eters £ and e. Based on a power spectrum generated 
using the modified version of CAMB, a TT data set was 
simulated by drawing from a ^-distribution with 11 + 1 
degrees of freedom. The data set and model are depicted 
in figure (J3j) , and the choice of values for £ = 0.0004 
and e = 0.01 are motivated by the Horava-Witten model 
[2lT ] of transplanckian physics In this model, unifi- 
cation occurs at the scale where a fifth dimension be- 
comes visible [TtJ , resulting in quite optimistic estimates 
of the energy scale where quantum gravity effects are 
manifested. However, the oscillations for this low value 
of £ corresponds to very small superimposed oscillations 
in the primordial power spectrum, as seen from figure 
[1] and [21 The remaining cosmological parameters used 
for creating the data set (eg £lb,T,n s ) were chosen as 
ACDM best-fit values. We continued by creating a sim- 
ulated TE, EE and BB data set using the same methods 
described here. 

The MCMC software package CosmoMC H2 was al- 
tered to include the parameters £ and e. Based on the 
likelihood for ideal noiseless simulated data described in 
[HI, the WMAP-likelihood code of CosmoMC was ad- 
justed to utilize the ideal data instead of the WMAP 
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spectrum used for estimating the transplanckian parame- 
ters. The results in might therefore have been partly 
due to the inaccuracy in the version of CMBFAST hjl 
used in this work. The accuracy of the numerical calcula- 
tions was boosted, which increased run time by a tenfold. 
The results can be seen in figure^ However, the accuracy 
boosting had little or no effect on the WMAP parameter 
estimation, as will be explained in the following section. 



FIG. 3: Simulated TT data based on £ = 0.0004 and e = 0.01. 



IV. RESULTS FROM WMAP3 DATA 



3-year data. The TT-only likelihood function is given by 

- 21og£(a; i ;p i ) = V(2*+ l)[ln &) + % - ll . (10) 
t L Li w J 

where CV are the theoretical TT power spectra and Ce 
TT data sets. The same method for creating and using 
ideal data with CosmoMC was employed in 24]. The 
errors in this ideal data set are only limited by cosmic 
variance, so no experiment can give any better estimate 
of the power spectrum. Thus, if transplanckian effects 
are not detectable in the simulated data, they cannot 
be detectable with any future high-precision data. We 
continued by employing the full simulated TT, TE, EE 
and BB data set using the likelihood presented in [25[ : 



-21n£ = 5^(2^+1) 
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For small values of the transplanckian values, the con- 
straints on both £ and e from the WMAP3 data turned 
out rather poor when running CosmoMC with all stan- 
dard cosmological parameters free. This can be seen in 
figure [5] and [71 consistent with the results in [l(| • For 
larger values of the transplanckian parameter £, the in- 
creasing amplitude of the primordial power spectrum ef- 
fectively constraints £ to the interval [0, 0.02] (recall from 
section HIl that for large £, the amplitude converges to 2). 
This can be seen in figure [S] This corresponds to a cutoff 
A - 0.002M p; ~ 10 16 GeV. Everything below this scale 
(£ > 0.02) should therefore be ruled out by WMAP3 
data. 



FIG. 5: Marginalized posterior for £ over a small interval 
surrounding the original input parameter £ = 0.0004 from 
the Horava-Witten model (red line). 



A. Increasing accuracy 
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FIG. 4: The ratio of modulated power spectrum with low and 
high accuracy, all 3 parameters boosted. 

Previous papers 0, HJ have pointed out that it is nec- 
essary to increase the accuracy of the calculated power 



FIG. 6: Marginalized posterior of £ over a larger interval. 
The red line corresponds to £ = 0.0004 from the Horava- 
Witten model. Notice how it is possible to estimate an upper 
bound: £ m ax ~ 0.02. Above this threshold, there are occa- 
sional "bumps" in the posterior where the oscillation reaches 
zero (sin ir = 0, see equation |9}. 

For physically reasonable values of £ and e (<C 1), the 
amplitude of the modulations are well below 1% of the 
total power spectrum. Also, the cosmic variance and the 



FIG. 7: Marginalized posterior for e. The red line corresponds 
to e = 0.01 from the Horava-Witten model. 




errors in the WMAP data are several orders of magni- 
tude greater than the amplitude of the modulations in 
the power spectrum, as can be seen from figure [8l This 
explains why any detection of transplanckian effects in 
the WMAP data fails for low values of £. 
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FIG. 8: The blue line shows the oscillations ACe in the 
primordial power spectrum for £ = 0.0004, 7 = 0.01 and 
e = 0.01. The red bars are the error bars in the WMAP data, 
while the green graph is the cosmic variance. Notice how the 
modulations in the low-£ power spectrum (blue line) are much 
smaller than either. The dotted blue line represents a larger 
£ ~ 0.01, which should be possible to detect, and gives upper 
constraints on £. 



V. RESULTS FROM IDEAL CMB DATA 

The WMAP3 data in CosmoMC was replaced by the 
optimal data set created from the power spectrum based 
on the Horava-Witten model. The ideal TT data are 
shown in figure [3] 



FIG. 9: Top left: Marginalized posterior for £ with mean 
value at 0.978 • 10 -3 . Top right: Marginalized posterior for 
£ with mean value at 0.515 • 10 -3 . Bottom left: Marginal- 
ized posterior for £ with no peak, almost uniform distribu- 
tion. Bottom right: Marginalized posterior for £ with several 
peaks. Notice how the peaks don't match, where the red line 
represents the Horava-Witten input parameter £ = 0.0004 for 
the simulated data. 

varying step length affect the end result: different initial 
parameters for the same parameter results in different 
distributions. This suggests that the likelihood function 
which decides the path of the random walkers traversing 
the cosmological parameter space has some troublesome 
properties. 




FIG. 10: Top left: Marginalized posterior for e with mean 
value at 0.12. Top right: Marginalized posterior for e with 
mean value at 0.037. Bottom left: Marginalized posterior for 
e with mean peak at 0.101. Notice how the peaks do not 
agree for different step lengths, where the red line represents 
the Horava-Witten input parameter e = 0.01 for the simulated 
data. 



A. Simulated TT data 

When performing parameter estimation of the trans- 
planckian parameters £ and e using the TT-data setup, 
the Monte Carlo Markov chain-method failed to converge 
to a stationary distribution. The resulting likelihoods 
for various £ and e are presented in figure [9] and 1101 re- 
spectively. Notice how different initial conditions and 



B. Simulated TT,TE,EE and BB data 

The marginalized posteriors can be seen in figure 1111 
Notice how r is restored due to the inclusion of simu- 
lated TE and EE data. The transplanckian parameters 
now reflect the ruggedness of the likelihood-function, and 
traces can also be seen in n s and A s . As in the TT-case, 
the original transplanckian input parameters are still not 
possible to pin down. 
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FIG. 11: Marginalized posteriors for the full TT, TE, EE and 
BB likelihood analysis using equation (fTTjl 



plots, so the best-fit value corresponds with the global 
minimum. The likelihood surfaces are riddled with local 
minima, effectively rendering the usual MCMC-method 
useless: the random walkers become "trapped" in the lo- 
cal holes, and will have problems converging to the cor- 
rect stationary distribution. Take especially note of the 
two-dimensional parameter space presented in figure 1161 
a landscape reminiscent that of an egg container. 
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FIG. 13: — log £ as a function of ^ for simulated transplanck- 
ian data. 
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FIG. 12: Marginalized posteriors (red lines) with smaller bins 
than in figure [U] for the full TT, TE, EE and BB likelihood 
analysis using equation (|llfl . The likelihoods are depicted 
with blue lines. Notice how the posteriors completely failed 
to converge to a distribution. 



Figure ITT1 and IT2l show how the MCMC process fails in 
converging, as the marginalized posterior show random 
walkers "trapped" in local minima. 



FIG. 14: — log £ as a function of £ for simulated non- 
transplanckian data. Notice how the non-transplanckian data 
is preferred to the simulated transplanckian data set. 
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VI. THE EXACT LIKELIHOOD FUNCTION 



FIG. 15: — log C between varying e and simulated data. 



In order to explain why the parameter estimation fails, 
we investigated the likelihood function in the case of 
noiseless data and different power spectra for fixed cos- 
mological parameters except the transplanckian £ and 
e. This was done by grid integration, and no MCMC- 
sampling was involved. The resulting likelihood sur- 
faces for continuous variations in £ and e are seen in 
figure H3] and [TS1 respectively. A figure of the exact two- 
dimensional likelihood surface for both £ and e is pre- 
sented in figure [ini Figures [T51 [T4"l [T5l and [TBI arc — log£- 



The shape of the exact likelihood is a result of the 
nature of the oscillations: equation ([9]) gives rise to 
a particularly nasty behaviour of the oscillations [29j |. 
The problem with the likelihood function is not easily 
solved. In principle, one could perform a smoothing 
of the likelihood, something which works well in a 
low-dimensional parameter space, but is hopelessly 
inefficient for a realistic number of parameters. 
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FIG. 16: — log C for varying both £ and e. The red ring 
represents the original input parameters £ = 0.0004, e = 0.01 
for the simulated data set. Notice how a random walker would 
be trapped in any of the other local minima. 




FIG. 17: -log Likelihood for varying £ and amplitude 



A. Decoupling the amplitude 

The amplitude and the frequency of the transplanckian 
oscillations are linked in the model we have considered 
so far. However, it is certainlypossible that they might 
vary independently, and in 0, [a, 0, Hi , the authors oper- 
ate with the amplitude of the transplanckian oscillations 
as a free amplitude. We repeated our analysis with fre- 
quency and amplitude decoupled. From figure [T"7l and [TBI 
it should be clear that freeing the amplitude does not 
really add any significance to the likelihood surface: the 
likelihood still gives the correct preferred amplitude. 



VII. BAYESIAN EVIDENCE 

From the logarithmic likelihood in figure [13] and I15[ it 
is obvious that it is impossible extract the values of the 
parameters of our model from the simulated data. How- 
ever, one might still ask the question of whether the data 
prefer the transplanckian model to the standard power- 
law primordial power spectrum. In a Bayesian framework 
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FIG. 18: -log Likelihood for varying £ and amplitude (differ- 
ent viewing angle) 

this question is answered by calculating the evidence ra- 
tio of the two models, see [26| for a nice introduction. 
Here we fix e = 0.01 and use the evidence ratio to see 
whether the data prefer the introduction of the additional 
parameter £. We assume that the transplanckian param- 
eter £ is not degenerate with any of the cosmological 
parameters, designated by the parameter vector 9. This 
might seem like a strong assumption, but we have checked 
that it is reasonable for the model and the parameters we 
consider in this paper. Also, the effects of primordial os- 
cillations have been investigated in [27j , concluding that 
degeneracies of oscillations with standard cosmological 
parameters are virtually non-existent. 

Let 9 denote the non-transplanckian cosmological pa- 
rameters, and 9 the best-fit KCDM parameters. The 
logarithmic evidence is given by 

A\nE = lnE -lnE Q (14) 

where E is the evidence of the non-transplanckian 
model. Then 

AhxE = {C}^) - (C) <e) (15) 

= (C(9, 0U() -C{6) (16) 

where tv(9, £) is the prior over [9, £) and ir(9) is the prior 
over 9. Thus, (C)x denotes the mean likelihood over 
the prior space X. The likelihood can be found in fig- 
ure [TH which is the exponential of the negative values 
in figure [131 We choose the prior to be £ £ [0,0.005]. 
This results in Ami? f» —1.5 which shows that the 
transplanckian model is slightly not preferred. Hence 
even perfect CMB-measurements over the whole sky 
does not give any significant evidence for transplanck- 
ian effects. When constraining the prior to reside in 
£ £ [0,0.0005], then A In E « 0.15 and the transplanck- 
ian model is more or less equally preferable to the non- 
transplanckian model. This is hardly surprising, as a 
continuous constraint of the prior Pr(£) = lim^^o[0,£] 
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FIG. 19: Likelihood £(£) used in calculating the significant 
evidence. 

results in A\nE = \nE° — \nE° = 0, which corresponds 
with the standard non-transplanckian model. The evi- 
dence never grows larger than ~ 0.15 for any prior con- 
straint on £. 



VIII. ALTERNATE MULTILEVEL SAMPLING 

There exists a few multilevel sampling methods be- 
sides the standard MCMC/Mctropolis-Hastings algo- 
rithm that can tame the ruffled likelihood through a se- 
ries of optimizations. Reference [28| contains a detailed 
introduction to these methods, most notably simulated 
annealing, umbrella sampling and the simulated temper- 
ing method. We proceed by a rough discussion on how 
some of these alternate sampling methods work when in- 
vestigating the 1-dimensional likelihood presented in fig- 
ure HU 



A. Metropolis-Hastings method, ID 



complicated likelihood surface prevents the posterior to 
converge to the correct shape. 

B. Modified Metropolis-Hastings method, ID 

Figure [3T] shows the posterior of 10 000 walkers after 
^20 000 steps each, using the Metropolis-Hastings algo- 
rithm. In addition, for each step there is a 0.01% chance 
for a given walker to perform a "large" jump to a new, 
(uniformly) random location. The walkers will therefore 
tend to accumulate at the correct location (where the 
-log likelihood is lowest, £ = 0.0004). This is a slow 
method, and needs a fair amount of steps to give results. 
It is also not especially "robust" , as the accumulation of 
walkers at the lowest point is dependent on fine-tuning 
the large-step rate (here, 0.01%). 




FIG. 21: The posterior of 10 000 MCMC walkers evaluated 
after 22 000 steps each in 1 dimension, with initial conditions 
£o = 0.0015. The posterior shows a clear contribution around 
the correct value of £ = 0.0004. 



C. Simulated annealing, ID 
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FIG. 20: The posterior of 10 000 MCMC walkers evaluated 
after 22 000 steps each in 1 dimension, with initial conditions 
£o = 0.0015. None of the walkers reached the correct value of 
£ = 0.0004. 

Figure [20] shows the posterior of 10 000 walkers after 
~ 22 000 steps each, using the Metropolis-Hastings algo- 
rithm. The likelihood is shown in figure 1131 In theory, 
the walkers should converge to a stationary distribution 
with the tallest peak located at £ = 0.0004, but due to 
the ruffled likelihood this will take forever. The initial 
values for the walkers was here £o = 0.0015. This reflects 
what we have seen in figure [H [TO1 [TT] and [TU : that the 
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FIG. 22: The result of 10 000 MCMC walkers evaluated after 
6 000 steps in 1 dimension, with initial conditions £o = 0.0015 
(not a posterior). The temperature was slowly decreased un- 
til reaching T = 0, where the correct value £ = 0.0004 was 
obtained. 

A description of the simulated annealing (SA) method 
can be found in [28|. As in the previous examples, 10 
000 random walkers were initialized at £o = 0.0015. The 
distribution spread out evenly as the SA temperature 
T was high, and started converging to different minima 
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as T was lowered. As T — > 0, a good portion of walkers 
ended up at £ = 0.0004, as seein in figure [22] This 
shows that the SA method can be successfully used to 
determine the correct minima for this type of likelihood 
landscape. 

However, the implementation and investigating of 
these alternative methods in CosmoMC is a major un- 
dertaking which we hope to pursue in future work. 

IX. CONCLUSION 

We have shown that standard MCMC methods are not 
suitable for detecting transplanckian oscillations in the 
CMB. The 3-year WMAP3 data alone gives weak upper 
constraints on the transplanckian parameters. With the 
simulated noiseless data, the results are more ambiguous. 
In theory, the estimation of the transplanckian parame- 
ters should be possible, but the messy likelihood func- 
tion renders conventional MCMC methods of little use. 
This is not due to the amplitude of the superimposed 
oscillations in the power spectrum alone, but rather re- 
flects the chaotic behaviour of the oscillations for varying 
transplanckian parameters £ and e. Small variations in 
the parameter space thus correspond to large variations 
in the likelihood function, as presented in figure [13] and 
M As in we conclude that CMB data are highly 
sensitive to these modulations in the primordial power 
spectrum, but it is virtually impossible to perform a pa- 
rameter estimation from WMAP3 data or high-precision 
simulated data using conventional MCMC methods. We 
have also shown, via Baycsian analysis, that even a high- 
precision measurement of the full CMB sky docs not give 
any significant evidence for transplanckian effects. It 
should however be possible to implement the SA-method 
in CosmoMC when using simulated high-precision data 



in order to detect the transplanckian parameters. 

We have in this paper shown that it is very unlikely 
that CMB-data could support the establishment of trans- 
planckian effects, even for models where quantum grav- 
ity is relevant on relatively low energy scales, like in the 
Horava-Witten model. There are several reasons for this 
lamentable conclusion: 

• A low quantum gravity scale, for example the value 
£ = 0.0004 chosen in our simulations, gives a very 
low amplitude of the oscillations in the primordial 
power spectrum. Hence, cosmic variance diffuses 
the (assumed noiseless) signal. 

• Evidence can maximally rise above 0.15, show- 
ing that the transplanckian model is only very 
weakly preferred when assuming a Horava-Witten 
universe. 

• Even large (and unphysical) values of £ will result 
in good fits for singular points due to the behaviour 
of the oscillations (they periodically reach zero 
when sin x = 0) 

• The conventional MCMC-methods in CosmoMC 
fails, and we have no current computational 
methods implemented for solving the problem of 
transversing the extremely complicated likelihood 
surface. 
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